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We present a new general mechanism to generate curvature perturbations after the end of the 
slow-roll phase of inflation. Our model is based on the simple assumption that the potential driving 
s ! ' inflation is characterized by an underlying global symmetry which is slightly broken. 
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I. INTRODUCTION 

\^ ' One of the most successful predictions of the inflationary theory, the current paradigm for understanding the 
evolution of the early universe |l| , is the redshifting of quantum fluctuations of the field driving inflation - the inflaton 
- beyond the Hubble radius, leading to an imprint on the background scalar (density) and tensor (gravitational waves) 
metric perturbations 0, IE EL IE IE EE IE IE \Es El E3 that subsequently seeds structure formation. 

For simplicity, most inflation models assume that there is only one scalar field involved in the dynamics of inflation. 
This is also the case when the mechanism of converting the energy driving inflation into radiation is considered. In 
this work we point out a qualitatively new effect that might arise if one relaxes the assumption of a single dynamical 
field. In a multi-field scenario in which the inflationary potential is characterized by a broken symmetry, the quantum 
fluctuations generated during the inflationary stage represent fluctuations in the initial conditions for the dynamics of 
the inflaton in the subsequent stage, thus implying that the background dynamics after the slow-roll phase has ended 
will differ in different regions of the universe. Since the background fields are coupled to the other fields into which 
& ' they decay, the fluctuations generated during the slow-roll phase will affect the subsequent decay process. 

The present work, assuming that the inflaton decay into other fields through the non-perturbative process of 
preheating |l3lll4| . is then aimed to understand whether isocurvature inflaton fluctuations, generated during the slow- 
roll stage, can lead to perturbations of the background metric through variations of the preheating efficiency. While 
the generation of curvature perturbations during the stages following the slow-roll phase has already been considered 
in some works HE O, E3, EH E3, E3, H3, , the present work is the first one to show that in a multi-field scenario 
a global broken symmetry of the potential is sufficient to yield curvature perturbations. Curvature perturbations 
produced through this mechanism can even represent the main source of perturbations to the background metric if 
the inflationary potential is such that the mass required to produce quantum fluctuations along the field trajectory is 
large, so that the latter result exponentially suppressed. 

The structure of the present work is the following. In Sec. ^ we obtain a general formula for the curvature 
perturbations generated from an inhomogeneous preheating efficiency related to the quantum fluctuations produced 
during inflation. Sec. IIIII presents an application of the general result obtained in Sec.|n]to the case of a broken U(l) 
symmetry. The conclusions are contained in Sec. IIVI 
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II. GENERAL RESULTS 

One of the main objectives in any particular preheating model is the calculation of the comoving number density 
of particles produced during the process, usually denoted by n x . In general, n x is a functional of the evolution and 

of the couplings of the preheat field x to the dynamically evolving field(s): n x = F[cj)(t)], where <j> — <j>\, 4>2 1 • • • , 4> n 
denotes the background inflaton fields that couple to x- Choosing a specific preheating model is then equivalent to 
specifying the functional F that relates <j)(t) to n x . 

In general, the dynamics of the background fields and of the scale factor is given by the solution of the system of 
coupled differential equations 

dV 

<p t + + — = 0, i = l,...,n, (la) 
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once the particular landscape of the potential V((f>) and a set of initial conditions {(p(to), <^(*o)} are specified. 

The main subject of the present work is to analyze how quantum fluctuations generated during the slow-roll stage 
of inflation may affect the efficiency of the subsequent preheating process. As inflation proceeds, quantum fluctuations 
will in fact cause the value of the background inflaton field to fluctuate in space about a mean value 4>o(t): 

$(t,x)=Mt)+5${t,x). (2) 

Denoting by to the epoch when inflation has ended but preheating has not yet commenced, it is possible to note 
that Eq. above shows that at each point in space the initial conditions 4>{to, %) that will determine the subsequent 
background field dynamics through Eqs. I|lallb[) are affected by the quantum fluctuations produced during inflation. 
Since the preheating efficiency is related to the dynamics of the background, it is then possible to conclude that 
quantum fluctuations produced during inflation may lead to fluctuations in the preheating efficiency through different 
background dynamics. 

Broken Symmetry. It is then necessary to point out that the mere presence of quantum fluctuations in the ini- 
tial conditions for the dynamics of the background during the preheating stage are not sufficient to yield different 
background evolutions leading to fluctuations of the preheating efficiency. If in fact the background potential V(<f>) 
is perfectly symmetric - that is V{<j>) = V(\<j>\) — then the fluctuations in the initial conditions will only lead to back- 
ground evolutions that are time translations of each other: in this case a simple rotation of the coordinate system in 
field space would again yield the well know case of a single scalar field. If the inflationary potential is characterized by 
a broken symmetry, on the other hand, then fluctuations in the initial conditions will lead to background trajectories 
that are not just time translations of one another. Two such background trajectories are shown in Fig. ^ f° r the 
case of a two-dimensional field space. Notice that the minimum distance to the origin in the trajectories are different. 
If the efficiency of the preheating process depends on the minimum distance obtained in the trajectory, then the 
preheating history will differ. 

Initial Conditions. It is important to distinguish here between the "initial conditions" for the background dynamics 
specified at the beginning of the slow-roll phase and those specified at time to, that is once the slow-roll phase has 
terminated and the preheating phase has not yet commenced. Since the present work deals with the particle production 
during preheating, the term "initial conditions" — and their fluctuations — refers here to the initial conditions for the 
preheating phase, that is to <fi(to,x). 

Considering the presence of the friction term 3H<fii in Eqs. (|laf> . it also seems reasonable to assume that the 
background dynamics during the preheating stage is mostly affected by the position that the background occupies 
in field space at the beginning of such a phase. Since the comoving number density of particles produced during 
preheating is a functional of the background trajectory, it is then possible to conclude that 

n x = n x [<p(t ), ...]. (3) 

Let's then turn to the preheating process and to the generation of curvature perturbations. As a first approximation, 
let's assume that preheating is complete and that the products of preheating are therefore the only particles populating 
the universe when the preheating stage has ended. 1 Neglecting the possible contributions stemming from non-adiabatic 



1 This assumption is quite important: since the preheat field \ is the only component present in the universe at the end of preheating 



3 



2 - 



2 

N 



2 0- 



c* 

■e. 



-3 



— 


t 

,•••'"*" 


1 B 1 


■ 
i 

•, — 




^ / /'''') 
I / /I 








i / / 1 


\ 






: II 


x i ■ 




— ; 


\ \ 1 

\ \ n 


/ 


| : — 




'■• If: 

\ (k 


0.5 "" 






\ "I k 


... i ** 






'HI '"• 


- 2 " /* 






/ /-•••... 








// 
i ' / i 


■■3 ' 

i i 


1 



-3-2-10 1 2 3 

0,/lO- 8 M p 

FIG. 1: Two background trajectories arising from fluctuations in the initial conditions in the case in which the background 
field space is two dimensional. Also shown are the equipotential contours (units of 10 _16 Mp). The potential used in this figure 
is the one given in Eq. 1161 and further analyzed in Sec. 11111 which is characterized by a global (7(1) broken symmetry. In this 
case, and for illustration purposes only, the symmetry breaking parameter has been arbitrarily set to x = 0.15. 



pressure perturbations present during the preheating stage, an estimate for the curvature perturbation £ can be 
obtained considering the number density perturbation, 



Px 



(4) 



where the spatially flat gauge has been assumed and the proportionality constant a depends on the redshifting of the 
particle produced. The above expression then allows to obtain an estimate of the curvature perturbations produced 
during the preheating stage induced by the fluctuations in the initial conditions <5</>(io) present at the beginning of such 
a stage because of the preceding inflationary stage. To proceed further it is then necessary to note that the coordinate 
system chosen to express the potential is not necessarily the one suited to analyze the perturbations arising during 
inflation, since in this coordinate system adiabatic and isocurvature perturbations are not decoupled. Adiabatic 
perturbations are most commonly thought to be the ones dominating the energy density during the inflationary 
stage, while the isocurvature ones are usually considered not to affect the energy density. Recalling the work of by 
Gordon et al. [2j|, it is also possible to note that given the inflationary trajectory (f(t), adiabatic perturbations 

correspond to perturbations along the direction tangent to the trajectory 4>{t), while isocurvature perturbations 

correspond to perturbations in the hyperplane orthogonal to 4>(t). This point can be intuitively understood once Eq. 
(|lafl is considered: since the motion of the background is driven by the gradient of the potential, it is reasonable 



this automatically ensures that only curvature perturbations are present at that point. On the other hand, if the preheating process 
doesn't turn all the energy initially stored in the </> field into the \ field, isocurvature perturbations can also result. 
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to expect that the trajectory will be tangent to the gradient. This in turn means that perturbations perpendicular 
to the trajectory will necessarily be isocurvature perturbations since they lie on equipotcntial hypcrsurfaces (thus 
perturbing the field values but not the energy density). Perturbations along the direction of the trajectory, on the 
other hand, being orthogonal to equipotcntial hypcrsurfaces, will necessarily affect the value of the energy density 
and therefore correspond to adiabatic perturbations. To decouple adiabatic and isocurvature perturbation it is then 
necessary to rotate the field space coordinate system so that one of the unit vectors lies along the tangent to the 
trajectory. Let's first of all define the unit vector u\\ which points in the direction of field space parallel to the (tangent 
to the) trajectory by 

, - &>) (5 ) 



l#o)| 



which allows to determine the component of the perturbation parallel to the trajectory simply by 8<j>\\ — 8<f>(to) ■ uii. 
Using the latter, it is then possible to define the unit vector u± which is orthogonal to the (tangent to the) trajectory 
by 

8<j>(t Q ) - 8(j)\\u\\ 

U_L = — =: , (6) 

\8(j)(t ) - 5(/>i\u\\\ 

which then leads to the component of the perturbation orthogonal to the trajectory 8cf>± = 8<f>(to) ■ u±. A new coor- 
dinate system in field space has then been defined, and the perturbation 6<p(to,x) has been decomposed accordingly: 



5<j){to, x) = 8<t>\\u\\ + 8<f>±u±. (7) 

In this new coordinate system 8(f>u represents an adiabatic perturbation while 8<p± represents an entropy perturbation. 
It is now quite intuitive to note that if the potential is characterized by a broken symmetry, then while adiabatic 
perturbations of the initial conditions 8((>u will lead to background trajectories that are differing just by a time 
translation, entropy perturbations will lead to trajectories that substantially differ from each other and that will 
therefore produce variations of the preheating efficiency. 

With this redefinition of the coordinate system of field space it is then possible to estimate the variation of the 
comoving density of particles produced during preheating due to fluctuations in the initial conditions generated during 
inflation. Note in fact that 

— = A 5 h + 52 5(f>1 - = V ln "x ' 5( t> (to, x), 8 
n x 0(p± 

but since perturbations of the initial conditions parallel to the field velocity will simply lead to background evolutions 
that are time translations of each other it is possible to conclude that dln(n x )/d<fi\\ = and that the mechanism under 
analysis is thus able to convert entropy perturbations into adiabatic ones. 2 Also, it is possible to envision models in 
which the inflationary dynamics is such that the adiabatic perturbations 8<j>\\ are exponentially suppressed. Combining 
Eqs. |@J and J5J) it is therefore possible to conclude that an estimate of the curvature perturbations produced by the 
inhomogeneous preheating efficiency connected to fluctuations in the background field dynamics originated during the 
inflationary phase is given by 

dln(n Y ) . , 

C«a— p^x, 9 

where it is interesting to note that while the 8<fij_ factor is determined during the inflationary stage, the 
a[d\n(n x )/dcf>±} factor is determined by the details assumed for the preheating process (and the associated par- 
ticle theory). The general conclusion that really seems worth stressing though is that on rather general grounds this 
model allows the conversion of entropy perturbations into curvature perturbations. 

It seems important to stress once more that the "initial conditions" that are considered in this work for the 
background field dynamics are the initial condition that result once inflation has ended (that is when a becomes 



2 This fact is obviously related to the assumption that the inflaton is supposed to completely decay into — and only into — the \ field. 



5 



negative). This is because what does affect the preheating efficiency is the evolution history of the background (j>(t) 
when it oscillates about the minimum of its potential. The value of 5<f>± appearing in Eq. @ should then be evaluated 
at the beginning of the oscillatory phase of the background, after the slow roll phase has ended. Considering a specific 
k— mode, it is then possible to note that the value of the amplitude of the quantum fluctuations S<p±(k) is frozen after 
the corresponding wavelength has exited the horizon and that therefore it can safely be evaluated when the wavelength 
A = 2ir/k crosses the horizon. 3 Recalling the general definition for the power spectrum of a generic quantity da 



k 3 

Vscr(k) = r-d&Tfc 
2-k z 



(10) 



it is possible to see that the power spectrum and the spectral index of the curvature perturbations obtained through 
this mechanism are given by 



r c (k) = 



n — 1 



dln(n x ) 



dlnV( d In Vs<j> ± 



d\nk 



d\i\k 



(11a) 
(lib) 



From these expressions it is also important to point out that while the power spectrum is affected by the specific 
nature of the preheating process, the spectral index is affected only by the characteristics of the potential in the region 
where quantum fluctuations are stretched to superhorizon scales (which are reflected in the power spectrum of d<f>±). 



III. APPLICATION TO THE BROKEN (7(1) CASE 

Let's apply the previous general results to the case in which the scalar field landscape is described by two degrees 
of freedom, <j>\ and 02- In this case, it is useful to express the potential in terms of a complex field <f>, 

4> = 4>i + ifo = W l9 . (12) 

If the potential V(<f>i, <f>2) is characterized by an exact U(l) global symmetry, then at the end of inflation the 
trajectory in field space will be in the radial direction. 4 In this case a simple rotation of the coordinate system in 
field space would yield again the well known case of a single scalar field, which then implies that fluctuations in the 
angular component - which in this case corresponds to the previous u± direction - of the initial conditions would not 
affect the background dynamics. Let's then investigate what are the consequences on the preheating process of an 
inflationary potential characterized by a slightly broken U(l) symmetry 



A. Assumptions and Basic Results 



Following the notation of the previous section, the initial conditions for the background field trajectory can be 
specified by [</>i(io), <h(to)] or by [|</>o|, #o] (where, since there is no possibility of confusion, the subscript here refers 
to the initial conditions) and their corresponding time derivatives. As was argued in the general case, the fluctuations 
in the initial field velocities can be neglected. Furthermore, recalling the presence of the 3H(j)i damping term in the 
background equations of motion it is possible to argue that after a first transient the trajectory in field space will be 
mostly along the radial direction. 5 It is therefore immediate to identify the new coordinate system for the field space 
as 

My ~ \4> \, (13a) 
u±_ ~ 9q. (13b) 



3 The calculation of the amplitude of such quantum fluctuation, along with its power spectrum and the resulting power spectrum and 
spectral index for the curvature perturbation, is presented in the next section for the specific case of a parabolic potential with a broken 
cylindrical symmetry. 

4 This is due to the fact that while the radial acceleration has a source term from the potential, if there is a U(l) symmetry then the 
angular component has only the damping term 3H 8 arising from the expansion of the universe. The presence of the damping term then 
causes any initial angular velocity 8o to decay away. 

5 "Mostly" because the symmetry breaking term contribute a small source term to the angular velocity. 
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In the present case the comoving number density of particles produced during preheating n x will therefore be 
a function of the initial conditions: n x — n x [\<po\, 6q, ...]. It is then possible to apply Eq. @ above to produce an 
estimate of the curvature perturbations produced during the preheating stage caused by the fluctuations in the angular 
direction 66q present at the beginning of such a stage, 



du 



(14) 



The power spectrum and the spectral index of the curvature perturbation thus obtained can also be estimated 
applying Eqs. I|llalllb|) : 



V c (k) 



n-1 = 



dln(n x ) 



d0 o 

d In d In Vse 



d\nk 



din k 



(15a) 
(15b) 



To proceed any further in the calculation it is necessary to specify the two details that so far have been left 
completely general. The first detail pertains the actual form of the inflatonary potential. Since the U(l) symmetry is 
assumed to be slightly broken, we assume that it takes the simple form 



(1 



(16) 



where x represents a measure of the symmetry breaking. 6 The origin of a nonvanishing value of x may be gravitational 
effects which can strongly violate global symmetries [24|. In such a case, x is likely to be given by (some power of) 
the ratio of the fundamental energy scale in the problem to the Planckian scale. 

The second detail that needs to be specified is the actual preheating model, thus nailing the specific functional that 
connects the background dynamics to the comoving number density of particles. In the present work the preheating 
model assumed is the instant preheating model of Felder et al. |14j. It is in fact not so unreasonable to suppose 
that the preheat field \ is coupled to some other fields into which it can decay. Furthermore, this choice for the 
preheating model is also characterized by some computational simplicity since in this case it is possible to express 
n x as a function of the initial conditions imposed on the background dynamics without having to resort to heavy 
numerical simulations. 

The instant preheating model assumes that the inflaton field 
standard (and simplest) preheating interaction, L^ x — — \g 2x 
field ip by the interaction C x ^ = h'tptpx- 7 Depending on the value of the coupling constants, the process (j) — > x ~~ * 
can be very efficient and turn the energy density initially stored in the background field into fermions in a single half 
oscillation of the inflaton about the minimum of its potential. 

Applying the results first obtained by Kofman et al. [Tflj . it is possible to compute the comoving number density 
of particles produced during the first pass of the background inflaton about the minimum of the potential. Given 
the interaction Lagrangian, it is important to note that if the inflaton trajectory doesn't exactly pass through the 
minimum of the potential (located at the origin of the coordinates in field space) but at a minimum distance \4>*\, 



is coupled to the preheat field \ through the 
X i an d that the field x 1S a ^ so coupled to a fermion 



then the preheat particles generated will be characterized b 
density of x particles produced in this case is then given by 



an effective mass m x = g\<f>*\. The comoving number 



y an er 

mm 



(g\k\) 



3/2 



■ exp 



irg\(f), 



21 



(17) 



where t„ denotes the instant in which the inflaton reaches the minimum of the potential V along its trajectory and 
and |</>*| respectively denote the field velocity and distance from the origin at such an instant. 



6 Recalling the fact that during preheating the value of n x is an adiabatic invariant and that it changes only when the background field 
is located in a small region of field space surrounding the the minimum of its potential, it is possible to note that Eq. 1 Kit represents 
quite a general choice since in such a region any potential can be well approximated in this form. On the other hand the spectrum 
of the initial condition perturbation S9o will depend on the form taken by the potential in the region where it drives inflation. The 
specification of Eq. 1161 above then can be considered general as far as the estimation of the d \n(n x ) / cLBq factor is concerned, but it is 
not general at all once the estimate of 89o(k) is considered. 

7 We suppose the fermion field ip to be massless, but this assumption is not crucial. 
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B. Estimate of the Curvature Perturbations 



Estimation of the curvature perturbations are not so complicated. Given the interaction Lagrangian connecting the 
preheat field x with the fermion field ip, the decay rate for the perturbative process x ~ * iprf is given by T x ^^ = 
h 2 m x /8ir = h 2 g\(j)\/8Tr. Since the decay rate for this process increases as the inflaton moves away from the minimum, it 
is not so unreasonable to assume that the whole process — * \ — > ip may be completed in a single half oscillation and 
that at the end only the fermions ip are going to be present. In the spatially flat gauge, the curvature perturbations 
will then be given by 



-H 



Sp^ 
Pi> ' 



(18) 



where p^ is the energy density of the fermion field. 

Recalling the form of the interaction Lagrangian that connects the inflaton and the preheat fields, it is also possible 
to note that when the inflaton moves away from the minimum of its potential it will endow the preheat field with an 
effective mass m x = g\<p\ which will quickly grow, rendering the preheat field nonrelativistic. Its energy density will 



therefore be given by p x 



i x g\<j>\. As a first approximation, let's suppose that all the x particles decay in 



one single instant t\ when the background inflaton field value is |0(<i)|. Then 



which in turn implies 



(19) 



(20) 



where the constant a depends on the redshift of the ip particles. In the present case, the ip particles are assumed to 
be massless so a — 1/4. From Eq. (|17fl . it is then straightforward to compute 



irg\(p*\ 2 \ <5|0*| 2irg\<p*\ 2 6\cp*\ 



\<P* 



|0»| 10* 



(21) 



To connect this expression with fluctuations of the initial conditions in the angular direction it is then necessary 
to express |0*| and |0*| as functions of the initial conditions. This is not such a complicated task given the fact that 
if we may neglect the expansion of the universe on the time scale of the single oscillation during which preheating 
occurs, we may exactly solve for the background dynamics (for more details, see Appendix 0. Expressing the initial 
conditions in polar coordinates as [|0o|,#o], the two parameters are approximately given by 



\M\M,0 o ;x)\ p 
|0*(|0o|,0o;z)| p 

As an intuitive check, it is possible to note that as x 



I0o i 

2V2 ' 



sin(26> 



x sin 



))■ 



(22a) 
(22b) 



and the symmetry is unbroken, Eqs. I|22al22b() yield the 
correct asymptotic behavior: |0*| vanishes because all background trajectories pass through the origin, and the #o 
dependence of |0*| also disappears. 

Given Eqs. I|22al22b|) . it is straightforward to compute their relative variations as a function of the variation in the 
initial condition of the angle S9q: 



= 2 



cos(26» ) 
sin(20 o ) 
x sin(2#o 



o = 2cot(20 o )*0o, 



2 [l-xsin 2 (6» )] 



(23a) 
(23b) 



Use of the above expressions in Eq. l(2T|) then yields 



5n, 



f(0 )56 a = 

sin(26> 



7Tffl0* 
10*1 



x sin(2#o 



47T3|0* 



2[l-xsin 2 (0 o )] 
3 Tr 3 \(po\gx 



■ cot(26» ) 



^l-xsm 2 (e ) l4^/l-xsin 2 (6» ) 



8to 



a;sin 2 (26» ) 
2[1 -xsm 2 (6» )] 



4cos(26» ) 



(24) 
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FIG. 2: Values of f(6o) obtained through numerical simulation and using the approximate expression given in Eq. 12411 for a 
value of the symmetry breaking parameter of x = 0.05 and for a coupling constant g — 0.01. 



It is interesting to note that while the second and third term in the curly brackets are coming from the exponential 
present in Eq. I|17|l , the first one comes from the multiplicative term. This means that for very small values of the 
symmetry-breaking parameter x the exponential suppression appearing in Eq. (|17() does not apply since the trajectories 
of the background field all pass extremely close to the origin of the coordinate system and the \ particles generated 
are almost massless. For larger values of the symmetry-breaking parameter the exponential suppression sets in and 
therefore the last two terms become crucial in determining Sn x /n x because small variation in the angle can lead to 
significant variations in the suppressing exponential. A plot comparing the expression obtained for f(9o) with the 
results of a numerical simulation is given in Fig. [21 The final expression for the curvature perturbations £ produced 
through this mechanism is thus given by 



x sin(2#o 



ir 3 \(j> \gx 




x sin 2 (26<o ) 



2[1 



x sin 



4cos(26> ) 



(25) 



where it is furthermore possible to note that in the x — ► limit £ vanishes. This fact is consistent with the point 
raised above regarding the case of a perfect symmetry: if the symmetry is unbroken then entropy perturbations will 
only lead to background evolutions that are time translations of one another and therefore the term d\n(n x )/ddo 
vanishes because n x doesn't depend on the angular initial condition. Since x plays a crucial role in determining the 
overall scale of the density perturbations, it would be interesting to investigate its value in realistic models, e.g. in 
those cases in which the breaking of the global symmetry is due to gravitational effects. 

Let's then proceed to compute the power spectrum and the spectral index of the curvature perturbations. It has 
already been argued that the amplitude |(5#o(fc)| can be evaluated at horizon exit. Letting Hk be the value of the 
Hubble parameter when the wavelength A = 27t/fc crossed the Hubble radius during inflation, it is possible to show 
that for the potential of Eq. (|16|l the square of the amplitude of the quantum fluctuations is given by 



\S9 (k)\ 2 = H * 



2fc 3 |0| 2 



aV" 1 

aH J 



1 — 2xrji cos 



'"''"l^f 



(26) 



where rji = m 2 /2>H 2 -C 1. Eq. (|26() shows that the symmetry breaking induces an very small correction to the ordinary 
flat power spectrum 



( H k 



V 2tt| 



aH J 



2;,i 



1 - 2riixcos z (9 a )ln 

1 aH 



(27) 



Now note that the factor f(9o) is uniquely determined by the initial conditions on the angle and by the specification 
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of the potential. The power spectrum of the curvature perturbations can then be obtained by 



V c (k) 



16 



Pse (k) 



f 2 (0o) ( H k 
16 \2tt\4> 



aH J 



2m 



k 

1 - 2r?ixcos ;!! (6 , o) In [ — — 
an 



f(0 o ) (Jh_ 

16 \2tt\4> 



(28) 



where the explicit form obtained for f(8o) in this case [Eq. 124(1 ] has not been entered for sake of brevity. Eq. 128|) 
then shows that the power spectrum of the curvature perturbations generated through this process is flat to a very 
good degree. This last aspect can be further stressed by the calculation of the spectral index, which gives 



1 dlnV < 



2r]ixcos 2 (9 ) 



1 - 2r ? ixcos 2 (6» )ln(fc/a J ff) 



2 [1 -a;cos 2 (6» )] 



3iJ 2 ' 



(29) 



thus showing that a very small tilt is induced by the symmetry breaking of the potential and the angle of the 
background trajectory. 



IV. DISCUSSION 



The analysis of this paper suggests that the production of curvature perturbations due to a broken symmetry of 
the inflationary potential and the resulting inhomogeneous efficiency of the preheating stage may be a rather common 
phenomenon. It is, in fact, appropriate to stress that while the magnitude of the curvature perturbations produced 
through this mechanism will depend on the details chosen for the specific model (potentials, coupling constants, 
interaction Lagrangians, preheating mechanism, and so forth), the mere fact that the inflationary potential has a 
minimum characterized by a broken symmetry is sufficient to guarantee the generation of curvature perturbations 
during the preheating phase. This is because perturbations in a direction orthogonal to the field trajectory yield 
evolutions of the background that are not just time translations of one other (as would be the case with perturbations 
along the direction of the trajectory) but that might differ substantially. Such different background evolutions then 
necessarily lead to different preheating efficiencies, thus resulting in perturbations in the comoving particle number 
and energy densities. 

As it has been shown in Sec. II I II choosing a specific preheating model allows one to quantify the magnitude of 
the curvature perturbations produced by this mechanism and to assess whether these may or may not represent a 
dominant component with respect to the adiabatic perturbation produced during the slow-roll phase by fluctuations 
along the radial direction. In the present context, the choice of the instant preheating model of Felder et al. |l4j | has 
been made because it seems plausible that the preheat field x ma y be coupled to some other fields. Moreover, the 
nature of such process allows one to obtain convenient analytic estimates of the comoving number density of particles 
produced that are not affected by the stochasticity, related to the build up of preheat particles, usually present in the 
standard preheating models 13]. Nonetheless, it seems important to stress that the main conclusions of this work do 
not depend on the choice of the preheating process, but only on the fact that n x depends on the background history, 
which in turn depends on the initial conditions. 

Finally, it seems important to note that if the preheat field is coupled to, and decays into, one single field, then 
the effect of a broken symmetry of the inflationary potential is to convert isocurvature perturbations into adiabatic 
perturbations |2^] . In this sense, the above model resembles in spirit the curvaton model of Lyth and Wands [l8| j but 
does not require the assumption of an external field. 
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APPENDIX A: APPROXIMATE VALUES OF \<t>»(\<j> \,9 ;x)\ AND \<j>»(\<j> \, 9 ; x) 



First consider Eq. I|22a|) . It is straightforward to note that given the inflationary potential, Eq. I|16[) , the exact solu- 
tion for the background dynamics can be computed once the expansion of the Universe is neglected. The background 
field dynamics is then given by 



4>i(t) = \4>o\ cos(#o) cos(mi), (Ala) 

rr, 



TTit 

Mt) = |0o|sin(0o)cos ( ^_ ) . (Alb) 
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Letting A = 7rx/4, it is possible to linearize the trajectory by 



4>i (t) « -|0o|cos(6> o )sin(A)A 
Mt) « |0o|sm(0 o )sm(A)(l-A), 



(A2a) 
(A2b) 



with A € [0, 1]. Then solving for the minimum distance from the origin yields A = sin 2 (#o), which corresponds to a 
minimum distance of 



10* 



(A3) 



where the approximation sin(A) w A has also been used. A plot comparing the analytical approximation obtained in 
this way with a numerical simulation is shown in Fig. |31 

Not let's then turn to Eq. I)22bl) . It is a basic fact of a simple harmonic oscillator that = m$ where $ is the 
amplitude of the oscillation and the mass m is given by the second derivative of the potential along the trajectory 



m 



d 2 V 



(A4) 



Here, / parametrizes the trajectory. To obtain a good estimate for \(j)*\ let's first assume that all trajectories pass 
through the origin. We define 



d 2 v 

d 2 v 

d<j>l 



l + x 



If the trajectory makes an angle (3 with the 4>% axis, it is then straightforward to show that 

sin 2 (/3)" 



mp = m 1 cos (/3) + m 2 sin (/3) = m 
which for small value of x reduces to 



cos 2 (/3) 



l+x 



mp w myl — x sin 2 (3. 



(A5a) 
(A5b) 

(A6) 
(A7) 
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tt/4 ti/2 3tt/4 u 

FIG. 4: Values of |0*|/|<^*| obtained through numerical simulation and using the approximate expression Eq. I|A70 for a value 
of the symmetry breaking parameter of x — 0.05. 



Assuming that the trajectory is perfectly radial, which is correct to a good approximation as long as the symmetry 
breaking parameter x is not too large, then (3 = 6q, and therefore 

|<£*| « 10*1^1 -xsin 2 O , (A8) 

where \<j> m \ — m\4>o\ depends only on the initial radial condition. 8 A plot comparing the analytical approximation 
given by Eq. I|A7() with numerical simulation is shown in Fig. 



APPENDIX B: CALCULATION OF THE POWER SPECTRUM AND THE SPECTRAL INDEX 

Start from the fact that 




where the two scalar fields are supposed to be uncorrelated. Then for a generic massive scalar field Uk we know that 
on superhorizon scales the amplitude of the quantum fluctuations is given by 

^M^V- <B3) 

where = 9/4 — m 2 a /H 2 and it is possible to define the parameter r\„ = m^/iH 2 ~ 3/2 — v a . In our case we have 
two fields, with masses — m and m^ 2 = m/Vl + x. We can then use Eq. I|B3|) to calculate \59 a {k)\ using Eq. 
(|B2|l . However since we're really interested in computing the power spectrum, we can go directly to its formula, which 



The fact that \<j>* \ = m|</>ol simply follows from applying conservation of energy to the background dynamics neglecting the expansion 
of the Universe term. 
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requires \S6o(k)\ 2 . In general, the power spectrum is defined by 

k :i 



Vsa^^\Sa k \ 2 , (B4) 



so that the first goal of our calculation becomes 



V se = J^\6e (k)\\ (B5) 



but using Eqs. IjBll IB2|I above we then have 



l^o(fc)| 2 = ^|#i| 2 + ^|# 2 | 2 . (B6) 
Now the only problem is that since the two fields have slightly different masses, we don't have \8<fii\ = \5(f>2\- Instead 



sm W ) 



k \ 2m 
— 1 +cos 2 (0 o ) 



aH I 



aH J 



27)2 



but 



(B7) 



>)2 



3H 2 (l + x) ZH 2 



(l-x) = ??i(l - x), 



so that 



\se a (k)\ 2 = H " 



2fc 3 | 



_^ 2 
aH 



But recall that x is very small, so that we can use the fact that a x w 1 + a; ln(a) which then yields 



l^o(fc)| 5 



HI 



1 — 2x77i cos 2 (^o ) In 



2fc 3 |(/>| 2 Vai?y 

We're now ready to insert this expression in the general definition of the power spectrum and then we get: 



(B8) 



(B9) 



(BIO) 



k 3 H 2 

Vs6 = 



2tt 2 2k 3 \(f> 2 \ \aH) 

H k y f k y 1 

27TI0IJ \aH) 



2'/i 



1 - 2i7ia;/3 2 In ( — ^ 
1 - 2r)ixl3 2 In 



~aH 



(Bll) 



We can then note that since cos 2 (#o) € [0, 1] and x -C 1 the spectrum we get is almost flat. Furthermore, the 
spectral index is given by 



71-1 



dlnP c 
d\nk 



2771 



2771 a; cos 2 ( 6>o ) 



1 - 277ia;cos 2 (6'o)ln(fc/ai/) 



2771 [1 -xcos 2 (6 )] = 2 [1 - xcos 2 (6»o)] 



(B12) 



where the denominator has been approximated to one since the 771X cos 2 ^) term should be very small. 
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